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We report on a search for the production of the Higgs boson decaying to two bottom quarks ac- 
companied by two additional quarks. The data sample used corresponds to an integrated luminosity 
of approximately 4fb~^ of pp collisions at y/s = 1.96 TeV recorded by the CDF II experiment. This 
search includes twice the integrated luminosity of the previous published result, uses analysis tech- 
niques to distinguish jets originating from light flavor quarks and those from gluon radiation, and 
adds sensitivity to a Higgs boson produced by vector boson fusion. We find no evidence of the Higgs 
boson and place limits on the Higgs boson production cross section for Higgs boson masses between 
lOOGeV/c^ and 150GeV/c^ at the 95% confidence level. For a Higgs boson mass of 120 GeV/c^ 
the observed (expected) hmit is 10.5 (20.0) times the predicted Standard Model cross section. 



I. INTRODUCTION 

The Higgs boson remains the only undiscovered parti- 
cle of the standard model (SM) of particle physics. It 
is the physical manifestation of the mechanism which 
provides mass to fundamental particles (T], Direct 
searches at the LEP collider have excluded a Higgs bo- 
son mass rriH < 114.4 GeV/c^ at 95% confidence level 
(CL) [3II , while the Tevatron collaborations have excluded 
a Higgs boson mass between 163 GeV/ and 166 GeV/c^ 
at 95% CL The Tevatron collaborations have re- 
ported a preliminary update which extends the exclu- 
sion region for a Higgs boson mass between 158 and 
173 GeV/ [5)]. Global fits to precision electroweak mea- 
surements set a one-sided 95% CL upper limit on mn at 
157GeV/c2 [61. 

This article presents the results of a search for the 
Higgs boson using an integrated luminosity of 4fb^^ of 
pp collision data at y/s — 1.96 TeV recorded by the Col- 
lider Detector at Fermilab (CDF H). We search for a 
Higgs boson decaying to a pair of bottom-quark jets (bb) 



accompanied by two additional quark jets (qq') for Higgs 
mass 100 < mn < 150GeV/c^. This search is most sen- 
sitive to a Higgs boson with low mass, mn < 135 GeV/c^, 
where the Higgs boson decay to bb is dominant [tI . The 
two production channels studied are associated produc- 
tion and vector boson fusion (VBF). The associated pro- 
duction channel is pp — VH qq' bb, where V is a W/Z 
vector boson, which decays to a pair of quarks. The 
hadronic branching fraction of V to qq' is ~ 70% Q. 
In the VBF channel, pp — ?> qq'H — > qq' bb, the incoming 
partons each radiate a vector boson and the two vector 
bosons fuse to form a Higgs boson. 

Low-mass Higgs boson searches at CDF have concen- 
trated on signatures that are a combination of jets, lep- 
tons and missing transverse energy which help to reduce 
the backgrounds but the signal yields are small [9l-[l]|. 
The hadronic modes used in this search exploit the larger 
branching fraction and thus have the largest signal yields 
among all the search channels at CDF. The major chal- 
lenge for this search is the modeling and suppression of 
the large background from QCD multijets (referred to as 
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QCD for brevity). 

A previous letter on the search for the Higgs boson 
in the ah-hadronic channel was published using an inte- 
grated luminosity of 2 fb~^ ^1^. This article has lowered 
the expected limit by a factor of two: a factor of « \f2 
from doubling the analyzed data and a factor of 1.4 from 
improvements to the analysis which are discussed in this 
article. 



II. THE TEVATRON AND THE CDF II 
DETECTOR 

The CDF II detector, designed to study -pp collisions, 
is both an azimuthally and forward-backward symmetric. 
It is described in detail in Refs. [Tsl - fTsI and references 
therein. CDF II uses a cylindrical coordinate system in 
which the z axis aligned along the proton beam direc- 
tion, Q is the polar angle relative to the z-axis and </) is 
the azimuthal angle relative to the x-axis. The pseudo- 
rapidity is defined as r\ = — ln(tan0/2). The transverse 
energy is Et = EsmO. Jets are defined by a cluster of 
energy in the calorimeter deposited inside a cone of ra- 
dius Ai? = ^/(A(/l)^ -I- = 0.4 as reconstructed by 
the JetClu algorithm [la]. Corrections are applied to 
the measured jet energy to account for detector calibra- 
tions, multiple interactions, underlying event and energy 
outside of the jet cone [l^ • 

The data for this search were collected by two mul- 
tijet triggers. The first 2.8 fb~^ used a trigger which 
selected at least four jet clusters with Et > 15 GeV for 
each jet and a total Et > 175 GeV. This trigger was used 
in the previous result [12]. The remaining 1.1 fb^^ were 
recorded with a new trigger which selected at least three 
jet clusters with Et > 20 GeV for each jet and a total 
Et > 130 GeV. The new trigger improved the accep- 
tance for a low-mass Higgs boson by 45% at Higgs mass 
of 100 GeV/c^ and by 20% at Higgs mass of 150 GeV/c^. 
The improvement was mainly due to lowering the total 
Et criteria in the new trigger. However the gain in the 
signal acceptance of the new trigger was diminished af- 
ter the event selection criteria which are described in the 
next section. 



background. A b jet is identified by its displaced vertex, 
as defined by the SecVtx algorithm ^1^], or by using the 
probability that the tracks within the jet are inconsistent 
with originating from the primary pp collision as defined 
by the JetProb algorithm 19]. The final four jets are 
labeled as 61, 62, 9i, 92 where b{q) are tagged (untagged) 
jets and E^ji'"' > 

The signal/background ratio is enhanced by dividing 
the data into two non-overlapping ^-tagging categories: 
SS when both jets are tagged by SecVtx, SJ when one 
jet is tagged by SecVtx and the other by JetProb. For 
a jet tagged by both algorithms, SecVtx takes prece- 
dence as it has a lower rate of misidentifying a light flavor 
jet as a & jet. T he p revious 2fb~^ search only included 
the SS category I12!] and the addition of the S J category 
increases the signal acceptance by 36%. Other 6-tagging 
combinations, such as both b jets selected by JetProb, 
were not considered in this search as the relative increase 
in the background is much larger than that for the signal. 

The data are divided into VH and VBF candidates 
defined by the invariant masses of the 6162 pair, mtb, 
and the (7192 pair, mqq. VH candidates have 75 < 
TUbb < 175GeV/c2 and 50 < < 120GeV/c2. VBF 
candidates have 75 < rribb < 175GeV/c^ and mqq > 
120GeV/c-^. The typical mbb dijet mass resolution is 
^ 18% [11]. These VH and VBF signal regions are illus- 
trated in Fig. [TJ We search for Higgs bosons produced 
via VH and VBF exclusively in the VH and VBF signal 
regions, respectively. The division of events is based on 
the different kinematics of the two processes. The VH 
channel has two mass resonances: mbb from the Higgs 
boson decay and mqq from the V decay. The VBF chan- 
nel shares the same mbb Higgs boson mass resonance but 
there is no accompanying resonance for mqq. The q jets 
in VBF tend to have a large 77 separation which results 
in larger values of mqq. The cut of mqq > 120GeV/c^ 
optimizes the VBF signal over background ratio. The 
acceptance for VH and VBF events varies from 2% to 
3% for 100GeV/c2 < mH < 150GeV/c2. As VJf and 
VBF candidates are also split by the two 6-tagging cate- 
gories, there are 4 independent samples (channels) which 
are studied: VH-SS; VH-SJ; VBF-SS; VBF-SJ. 



III. EVENT SELECTION 

Events with isolated leptons or missing transverse en- 
ergy significance [Ts*] > 6 are removed to avoid any over- 
lap with other low-mass Higgs analyses at CDF II. The 
data are refined further by selecting events with four or 
five jets where each jet has Et > 15 GeV and \r]\ < 2.4. 
The selected jets are ordered by descending jet-ET and 
any fifth jet plays no further role. The scalar sum of the 
four leading jets' Et is required to be > 220 GeV, and 
exactly two of the four leading jets are required to be 
identified ("tagged") as bottom-quark jets (6 jet). The 
scalar sum Et cut reduces the contribution of the QCD 



IV. SIGNAL AND BACKGROUND SAMPLES 

The data are compared to a model of the signal 
and background composed of QCD, tt, Z{-^ bb/cc) 
+ Jets (Z-l-jets), single-top, W + bb/cc (M^-|-HF), and 
WW/WZ/ZZ (diboson) events. The signal and non- 
QCD backgrounds are modeled by Monte Carlo (MC) 
simulation. The VH and VBF production are gener- 
ated by PYTHIA [13], combined with a GEANT-based (2l| 
simulation of the CDF II detect or y2/M . The non-QCD 
MC is described in detail in Ref. [ij] and normalized to 
next-to-leading order cross sections. All the MC sam- 
ples include the trigger simulation and their trig ger ef- 
ficiencies are corrected as described in Ref. 113]. The 



5 



m^^ (GeV/c') 
200 --■ 
175 --■ 



75 
50 



(3 



TAG 



CONTROL 



VH Signal 
Region 



CONTROL 



■/'////// 



VBF Signal 
Region 



m^^ (GeV/c') 

FIG. 1: mf,h — m^^ plane: This plane illustrates the VH 

and VBF signal regions used to select VH and VBF 
candidates. The Tag region is used to derive a model of 
the QCD background. The Control region measures 
the systematic uncertainty of the QCD model. 



tagged events) to predict the distribution of events with 
exactly two 6-tagged jets (double-tagged events). For 
each single-tagged event, the TRF gives the probability 
of each additional jet, called a probe jet, to be a second 
6-tagged jet. The TRF is parameterized as a function of 
three variables: the Et and ?/ of the probe jet and AR 
between the tagged b jet and probe jet. The choice of 
variables used to parameterize the TRF is motivated by 
the kinematics of the QCD background and the charac- 
teristics of the 6-tagging algorithms. As the behavior of 
the SecVtx and JetProb 6-tagging algorithms are not 
identical, there is a TRF for SS and another TRF for 
SJ. The TRF is measured using jets in the Tag region 
(Fig. [U, defined as m,, < 45GeV/c^, nibb < 50GeV/c^ 
and nibb > 200GeV/c2 which is not in the VH and VBF 
signal regions. 



VI. JET MOMENT 



TABLE I: Expected number of non-QCD background 
and VH/VBF signal with observed number of events for 
the four channels. Statistical and systematic 
uncertainties are combined in quadrature where 
systematic uncertainties dominate. The number of 
VH(VBF) events are exclusive to the VH(VBF) 
channels. The difference between data and non-QCD 
are assumed to be QCD. 





VH-SS 


VH-SJ 


VBF-SS 


VBF-SJ 


tt 


281.7 ± 45.6 


115.3 ± 19.9 


177.3 ± 28.7 


75.7 ± 13.1 


Single-top 


44.1 ± 7.1 


17.7 ± 3.1 


17.2 ± 2.8 


10.0 ± 1.7 


Z + Jets 


127.5 ± 65.8 


55.4 ± 28.8 


135.0 ± 69.7 


62.9 ± 32.7 


W + HF 


27.9 ± 14.4 


12.0 ± 6.2 


4.8 ± 2.5 


3.3 ± 1.7 


Diboson 


11.4 ± 1.6 


8.5 ± 1.3 


5.3 ± 0.7 


3.8 ± 0.6 


Total 


492.6 ± 81.7 


208.9 ± 35.7 


339.6 ± 75.5 


155.7 ± 35.3 


Non-QCD 










Higgs Signal (mj, = 120GcV/c'') 






VH 


7.8 ± 1.0 


2.9 ± 0.4 






VBF 






3.2 ± 0.4 


1.2 ± 0.2 


Data 


16857 


9341 


17776 


9518 



QCD background shape is modeled by a data-driven tech- 
nique developed in Ref. [l^l and described in detail be- 
low. The expected signal yields of the four channels 
are 7.8( VH-SS), 2.9( VH-SJ), 3.2(VBF-SS), and 1.2(VBF- 
SJ) for niH = 120GeV/c^. The total backgrounds are 
about 17000( VH-SS), 9300( VH-SJ), 18000(VBF-SS), and 
9500(VBF-SJ). The background composition is -98% 
QCD (Table IJ). 



V. QCD MODELING 

The shape of the dominant QCD background is mod- 
eled using a data-driven method known as the tag rate 
function (TRF) and is described in detail in [l^. The 
TRF is applied to a QCD dominated data sample of 
events with at least one SecVtx 6-tagged jet (single- 



The VH and VBF (/-jets are mostly quark jets while 
QCD g-jets are a mixture of gluon and quark jets. As 
gluon jets, on average, tend to be broader than quark 
jets, any variable related to the jet width is an additional 
tool to discriminate the Higgs signal from QCD. 

In this article, we use the jet 4>{r]) moment, {(j)) 
iiv)) of 9 jsts which measures the jet width along 
the <t>(r]) axis. The jet (/) and rj moments are defined by 



(0) = 



E 



^towcr 



E; 



,jct 



(^A0(0towor 



0jot)) 



{v) = 



\ 



E 

towers 



^towcr 



E; 



ijct 



f _ V 

1 ^towcr ^jct I 



(la) 



(lb) 



where the (j) and rj jet moments are summed over the 
calorimeter towers forming the jet and depends on the 
tower-Er {E^^^"'), the ^^et-Er (£^^^*), the tower's </>(r/) 
position, (/)towor (?7towcr), and the jet's 0(77) position, 
0jct(?yjct)- The function A(^(0towor, 0jet) in Eq. [Ta]is the 
smallest angular difference between </itower and 0jot. The 
jet moment is a measure of the jet's width. 

We checked whether the MC simulation of the quark 
jet moment matches the data. Gluon jet moments were 
not checked as the Higgs g-jet, modelled by MC, are 
mostly quark jets whereas gluon jets only appear in QCD 
which is derived from data. The hadronic W decay from 
tt — )■ bWbW blv + bqq' , where I is an electron or a 
muon, provides a source of quark jets. The event selec- 
tion from Ref. [25| was used to extract a tt data sample 
which is 86% tt. The complete sample composition is de- 
scribed in Ref. [l^l . The leading untagged jet pair whose 
invariant mass is 80 ± 30GeV/c^ is assumed to be the 
quark jets from the hadronic W boson decay. The same 
event and leading untagged jet pair selection is applied 
to ti MC to compare with data. 
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The jet moment depends not only on the parton initi- 
ating the jet but also on the -E^*, ?7jetj and the number of 
primary vertices in the event (A^tx) which are not guar- 
anteed to be the same for data and MC. The dependen- 
cies are removed by rescaling the measured jet moment 
to a common reference of Ei^^ —50 GcY / , '7jct=0 and 
iVvtx=l, as measured in data. The rescaling for (0) is 
performed using 

(0)r)ata = {4')Data X 

fLtaiE'r = 50 GeV/c^ 77j,t = 0, ^vtx ^ 1) ^^^^ 

/Data(^T'''7jct,^Vtx) 

(0)mc = {4')mc X 

fLtajE'r' = 50 GeV/c^ 77jet ^ 0, A^vtx ^ 1) ^^b) 

/mc(^t*' '7jot,^Vtx) 

where foatai^T^ ^Vjct, Nytx) and ff.j^{E^^\r]jct, N^t^) 
are the (0) parameterizations for data and MC, re- 
spectively. {4>) Data{{4>) Mc) are the measured ((/)) for 
data(MC), and {4')'Data{{'t')'Mc) rescaled values. 

(77) are rescaled in a similar way but has a separate {rf) 
parameterization for data and MC. 

After the measured jet moments are rescaled, the MC 
required an additional shift in (j) and 77 of ~ -\-2% to agree 
with the data. Half of this offset was used as an esti- 
mate of the systematic uncertainty of the MC jet mo- 
ment. Figure [3] compares the jet moments of data to the 
simulated ti signal and background MC which are in the 
same fractions as measured in data. Only after applying 
all corrections does the MC agree with the data. 

As an additional check, the {4>)'mc ^^'^ {''1)'mc 
MC was compared with VH and VBF MC. The aver- 
age jet moments of the MC samples were expected to be 
identical to the g-jets from tt as the Higgs signal are just 
quark jets. The jet moments from the VH sample agreed 
with tt. However there was a disagreement of 5% be- 
tween VBF MC and tt MC for jets with hi > 1.1. Half 
of this difference was used as an additional systematic 
uncertainty for the VBF jet moment. 

VII. NEURAL NETWORK 

The large background precludes the use of simple vari- 
ables, such as rubb, to search for a Higgs boson signal. 
An artificial neural network (NN), from the TMVA pack- 
age 25], is trained to search exclusively for VTf(VBF) 
Higgs bosons in the VH(VBF) signal region. The NN 
is trained using VH(VBF) signal and TRF QCD predic- 
tion as background. As the kinematics for VH and VBF 
Higgs signals are different, a dedicated NN for each sig- 
nal is trained. The NN training variables for the VH NN 
are iribb, rugq, rescaled qi 4> moment rescaled qi 7] 

moment (?7i)', rescaled (72 <!> moment {4>2)\ rescaled q2 p 
moment (772)', the cosine of the helicity angle cos 9* [26| . 



the cosine of the leading jet scattering angle in the four 
jet rest frame, cos ^3 ,27.] , and x which is a measure 
of whether both the 6-jet pair and g-jet pair are from 
a Higgs boson and V decay, respectively, x is defined 
as the minimum of xw and xz where xw is defined as 
Xw — \/{Mw — MqqY -\- [Mh — MbbY and a similar ex- 
pression exists for xz- For the VBF channel, the neural 
net inputs are rubb, mqq, {<j)i)' , {r]i)', (02)', and (772)'. 

The two 6-tagging categories have similar kinematic 
distributions which allows the same NN to be used for 
SS and SJ events. The NN is trained with SS events as 
it has the better signal/background ratio. 

Before training the NN, the TRF QCD modeling was 
verified by comparing the shapes of the NN training vari- 
ables constructed from single-tagged events, after apply- 
ing the TRF, and double tagged events from the TAG 
region. The TRF was able to reproduce the shapes of 
all the NN training variables except niqq and the jet mo- 
ments. The TRF was corrected using a correction func- 
tion for each mismodeled variable. The correction func- 
tion was constructed from the fitted ratio of the observed 
double tagged shape in the Tag region to the TRF pre- 
diction in the Tag region. The largest correction value 
was 2%. Figures |3] and ID show distributions of the NN 
training variables of VH and VBF signal, corrected TRF 
QCD and double tagged data for the SS 6-tagging cate- 
gory. The corrected TRF QCD follows the shape of the 
data for all variables. The TRF predictions for SJ were 
validated in the same way. 

We search for a Higgs boson of mass 100 < mn < 
150GeV/c^ at 5GeV/c^ intervals. As rubb is one 
of the NN training variables, which varies with differ- 
ent Higgs mass hypotheses, the Higgs search sensitiv- 
ity can be improved by training the NN at different 
Higgs masses. There is a separate yH(VBF) NN trained 
at niH =100GeV/c2, 120GcV/c2, and MOGeV/c^. 
For Higgs mass hypotheses between lOOGeV/c^ and 
llOGeV/c^, the NN trained with = lOOGeV/c^ is 
used. Similarly, we use the nin — 120GeV/c^ trained 
NN to search for a Higgs boson between 115GeV/c^ 
and 130GeV/c2 and the ttih = 140GeV/c2 trained NN 
to search for a Higgs boson between 135GeV/c^ and 
150GeV/c2. 

Figure [5] shows the NN distributions for VH and VBF 
for a Higgs mass of 120GeV/c^. The NN returns a more 
negative (positive) score for background (signal) events. 
As the QCD background is large, QCD subtracted NN 
distributions are also shown. 



VIII. SYSTEMATIC UNCERTAINTIES 

We estimate the effect of systematic uncertainties by 
propagating uncertainties on the NN input variables to 
the NN output. We consider both variations on the nor- 
malization and shape of the NN output. The system- 
atic uncertainties which affect the normalization of the 
Higgs signal and non-QCD backgrounds are: jet energy 



7 



350 



c 

UJ 300- 
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FIG. 2: Comparisons of the jet moments of tt data to tt and W^+HF, Z+jets, single top, diboson, misidcntified 
6-jets (M^+HF and others) MC. Only after applying all corrections does the MC agree with data. 



scale (JES) [T^l, parton distribution function (PDF) , b- 
tagging scale factor between MC and data, initial and 
final state radiation (ISR/FSR), trigger efficiency, inte- 
grated luminosity and cross sections [5| . The Higgs signal 
cross section uncertainty is ~ 5% ^ but is not included 
as its effect is negligible. 

The uncertainties which affect the Higgs signal NN 
output shape are the JES, ISR/FSR, and jet moments. 
Sec. IVII defined the jet moment uncertainty as half of the 
offset required to correct the MC. The Higgs signal jet 
moment NN shape uncertainty is defined as the differ- 
ence of the Higgs signal NN shape using the nominal jet 
moment correction and the Higgs signal NN shape using 
half of the jet moment correction. 

For the TRF QCD prediction, we consider two shape 
uncertainties: an interpolation uncertainty and correc- 
tion function uncertainty for m^q, ((/)), and {rj) variables. 
The interpolation uncertainty accounts for possible dif- 
ference in the TRF between the regions where it was 
measured (Tag) and applied (Signal) (Fig.[T|). An alter- 
native TRF was measured using events in the Control 
region, as indicated in Fig. [1] which is still background- 
dominated. The interpolation uncertainty is defined as 
the difference of the QCD NN shapes using the nominal 
Tag TRF and Control TRF. The correction function 
uncertainty for mqq is evaluated by deriving an alterna- 
tive correction function for nigq using events from the 
CONTROL region. This alternative nigq correction func- 
tion is applied to the TRF, instead of the nominal mqq 
correction function, and propagated through the NN. 
The difference in the QCD NN shape between using the 
nominal correction function and the Control region de- 
rived function defines the correction function shape un- 
certainty. The systematic uncertainty for (^), and (77) 
correction functions is evaluated in the same way. The 
QCD NN output varied at most by ~ 2 — 3% for each 
QCD shape systematic. The uncertainties are summa- 



TABLE II: Summary of systematic uncertainties. The 
largest change is quoted for the sources which have a 
shape uncertainty. 



Source Uncertainty 



Higgs and non-QCD uncertainties 




Integrated Luminosity 


±6% 


Trigger efficiency 


±4% 


PDF 


±2% 


JES 


±7% and shape 


b-tagging scale factor 


±7.6% for SS 




±9.7% for SJ 


ISR/FSR 


±2% and shape for VH 




±3% and shape for VBF 


(0)' and (77) ' 


shape for VH and VBF (< 10%) 


tt and single-top cross section 


±10% 


Diboson cross section 


±6% 


VK+HF & Z+Jets cross section 


±50% 


QCD uncertainties 




interpolation 


shape (< 3%) 


rriqq correction function 


shape (< 1%) 


(0) correction function 


shape (< 2%) 


(77) correction function 


shape (< 2%) 



rized in Table HIl 



IX. RESULTS 

The NN output distribution of data are compared to 
the background and we find no excess of events over the 
expected background. We calculate upper limits on the 
excluded Higgs boson cross section at the 95% CL for 
Higgs boson mass hypotheses 100 < niH < 150GeV/c^ 
at 5GeV/c^ intervals. The limits are calculated using a 
Bayesian likelihood method with a flat prior for the signal 
cross section. We integrate over Gaussian priors for the 
systematic uncertainties and incorporate correlated rate 
and shape uncertainties as well as uncorrelated bin-by- 
bin statistical uncertainties 1. The QCD normalization 
is a free parameter that is fit to the data. 
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FIG. 3: Distribution of variables used to train the VH NN. The signal consists of VH {mH = 120GeV/c^) SS events 
and the background consists of TRF predicted QCD SS events which have passed the VH candidate selection. All 
plots are normalised to unit area to compare shapes. After correction functions have been applied to mqq and the 
jet-moments, the TRF correctly predicts the shape of the double-tagged SS data for all variables 
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(a) (b) (c) 




FIG. 4: Distribution of variables used to train the VBF NN. The signal consists of VBF {mn = 120GeV/c2) SS 
events and the background consists of TRF predicted QCD SS events and which have passed the VBF candidate 
selection. All plots are normalised to unit area to compare shapes. After correction functions have been applied to 
mqq and the jet-moments, the TRF correctly predicts the shape of the double-tagged SS data for all variables. 



10 




<2 150f — I ^ ^ 

C - Data-QCD ^ QCD Uncertainty | | Z+jets 




-100 



-1.5 



-0.5 



0.5 1 1.5 
Neural Network Output 



(c) 



g> 150 
> 

lU 



100 



SO 



-50 



-100 



- Dala-QCD QCD Uncertainty |^ Z+jels 
I tf W+HF/Single Top/Diboson VBF x IOC 



-1.5 




-0.5 



0.5 1 1.5 2 
Neural Network Output 



(d) 



FIG. 5: NN distribution for VH-SS (a) and VBF-SS (b) for mn = 120GeV/c2. As the QCD background is large, 
data-QCD versions for yfZ-SS and VBF-SS are shown in (c) and (d), respectively. The VH and VBF distributions 
are scaled by a factor of 1000 and 100 for (a), (b) and (c), (d), respectively. 



The median of the 95% CL obtained from 10000 simu- 
lated experiments is taken at the expected 95% CL. The 
±lf7 (where a denotes the standard deviation) and ±2cr 
expected limits are derived from the 16th, 84th, 2nd and 
98th percentiles of the distribution, respectively. 

For niH = 120GeV/c^, the observed (expected) limit, 
normalized to the SM cross section, for the individual 
analysis channels are 11.9(25.6) for VH-SS, 43.4(51.8) for 
VH-SJ, 47.0(49.4) for VBF-SS, 93.7(132.3) for VBF-SJ, 
and 10.5(20.0) for the combination of these four channels. 
The combined channel limits for Higgs boson masses in 
the range between 100 - 150GeV/c^ are shown in Fig. [B] 
and summarized in Table Hill 

The observed limits for the individual search channels 



agree within la of their expected limit except for the 
VH-SS channel where we see a 2a discrepancy. The ob- 
served data in the VH-SS channel have a deficit in the 
high signal region of the NN. Since the VH-SS channel is 
the most sensitive, it has the strongest influence on the 
combined limit; thus the deviation of the observed limit 
from the expected limit in the VH-SS channel is similar 
to that of the combined limit. Figure [7] shows the ratio of 
the data to the expected background for the four analysis 
channels for the NN trained on a 120GeV/c'^ Higgs bo- 
son. All four channels show a ratio w 1 over the whole NN 
output range but the VH-SS channel has several points 
with a ratio of « 0.9 at the NN output of ~ 0.5; the 
most sensitive region of the NN output where the Higgs 
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TABLE III: Expected and observed 95% CL upper 
limits for the combined VH and VBF channels. The 
hmits are normahzed to the expected Higgs cross 
section. 



Higgs mass 


-2cr 


— 1(7 Median 


+l<y 


+2<y 


Observed 


(GcV/c^) 













100 


9.1 


12.8 


18.8 


27.2 


38.5 


10.1 


105 


8.7 


12.1 


17.4 


25.2 


35.8 


9.9 


110 


8.0 


11.7 


17.1 


24.5 


34.2 


10.2 


115 


8.8 


12.2 


17.8 


25.9 


36.9 


9.1 


120 


9.3 


13.7 


20.0 


28.5 


39.5 


10.5 


125 


13.5 


18.7 


27.3 


39.8 


57.0 


13.8 


130 


17.0 


24.4 


36.1 


52.8 


75.4 


17.2 


135 


19.6 


28.6 


41.9 


59.7 


82.7 


22.7 


140 


26.7 


40.7 


60.4 


86.6 


120.2 


35.2 


145 


43.4 


63.5 


95.7 


142.1 


205.3 


55.8 


150 


73.8 


109.9 


164.1 


240.3 


341.9 


101.0 



□ 



S 
o 



















— Observed 






-- Expected 






□ ±1a 






□ ±2a 









100 105 110 115 120 125 130 135 140 145 150 

Higgs Mass (GeV/c^) 

FIG. 6: Expected (dashed) and observed (solid) 95% 
CL normalized to the SM cross section for the combined 
VH and VBF channel. The dark (hght) band represents 
the la{2a) expected limit range. 



signal peaks. If the background was mismodeled, either 
the TRF has incorrectly predicted the QCD background 
or the NN was at fault. The VH-SS and VBF-SS chan- 
nels share the same TRF and the VBF-SS observed limit 
agrees with its expected limit. The VH-SS and VfZ-SJ 
channel share the same trained NN and the observed and 
expected limits in the WJ-SJ channel agree. Since nei- 
ther the NN nor TRF showed any evidence of corrupting 
the background prediction, it suggested the low ratio for 
VH-SS was likely to be a statistical fluctuation rather 
than evidence of background mismodeling. 



X. SUMMARY 

In summary, a search for the Higgs boson was per- 
formed in the all-hadronic final state and set observed 
(expected) limits of 10.5 (20.0) times the predicted Stan- 
dard Model cross section at 95% CL for 120GeV/c^ 
Higgs boson. The measurements presented in this ar- 
ticle has shown a factor of two improvement over the 
previous 2fb~^ result for the all-hadronic Higgs boson 
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FIG. 7: Ratios of the data to background for VTf-SS 
(a), VH-SJ (b), VBF-SS (c), and VBF-SJ (d) for the 
NN trained on 120GeV/c2 Higgs boson MC. The error 
bars of the data to background ratio are the statistical 

errors. The solid gray band is the ratio of the 
background systematic uncertainty to the background. 
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search [l^ . This article extended the 2^""*^ analysis by 
including the VBF channel, adding an additional algo- 
rithm to identify bottom-quark jets, adding an artificial 
neural network to separate signal from background which 
includes (</>) and (77) to distinguish gluon jets from quark 
jets, and by doubling the analyzed data set. CDF II con- 
tinues to collect more data and further improvements to 
the analysis technique will extend the sensitivity of the 
all-hadronic Higgs boson search. 
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